AD- A 185  023 


BTC  FILE  CORy 

***L-R-1048 


DEPARTMENT  OF  DEFENCE 

DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
MATERIALS  RESEARCH  LABORATORIES 
MELBOURNE,  VICTORIA 

REPORT 


DEPARTMENT  OF  DEFENCE 
MATERIALS  RESEARCH  LABORATORIES 


REPORT 


* 

5  MRL— R— 1048 


a 

a 

H 


I 


2D  MAGNETOSTATIC  CALCULATIONS  AND  IMPLICATIONS 
FOR  FERROMAGNETIC  MINE  CASES 


t 

P.N.  Johnston 


ABSTRACT 


: 


■3 

r 


I 


in  order  to  assess  the  effect  of  a  magnetic  mine  case  on  a  magnetic 
sensor  In  or  near  the  case,  some  2D  magnetostatic  calculations  have  been 
performed  using  the  boundary  element  method.  For  the  case  of  a  cylindrical 
tube  the  analytic  solution  to  the  problem  Is  presented  allowing  the  analytic 
and  numerical  solutions  to  be  compared.  Conclusions  are  drawn  about  the  ► 
likely  effect  of  a  ferromagnetic  case  on  a  sensor  In  or  near  the  case. 
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2D  MAGNETOSTATIC  CALCULATIONS  AND  IMPLICATIONS 


FOR  FERROMAGNETIC  MINE  CASES 


1 .  INTRODUCTION 

“  t' 

i 

The  presence  of  a  ferromagnetic  object  in  an  otherwise  uniform 
magnetic  field  causes  an  anomaly  (local  distortion  of  the  field).  The  anomaly 
can  be  used  to  signal  the  presence  of  a  ferromagnetic  object  e.g.  a  buried 
piece  of  ordnance.  A  time-varying  anomaly  can  be  used  to  signal  the  arrival  of 
a  vehicle  e.g.  sensors  in  the  road  connected  to  traffic  lights. 

The  presence  of  magnetic  materials  in  the  vicinity  of  a  magnetic 
sensor  will  alter  its  apparent  sensitivity  (its  real  sensitivity  is  of  course 
constant,  but  the  magnetic  field  distortion  makes  it  appear  changed)  to 
magnetic  field  changes.  This  poses  a  problem  when  considering  a  sensor  in  a 
mine  with  a  ferromagnetic  case.  Many  currently  available  mines  have  this 
arrangement . 


The  magnetic  sensor  of  ground  mines  (sea  mines  laid  on  the  sea 
bottom)  is  usually  the  sensor  used  to  determine  (i)  whether  the  mine  is 
actuated  or  the  ship-count  is  progressed  and  (ii)  the  position  along  the 
ship's  path  relative  to  the  mine  at  which  actuation  occurs.  Therefore  the 
apparent  sensitivity  of  the  magnetic  sensor  affects  the  accuracy  of  target 
ranging.  If  the  sensitivity  of  the  mine  is  variable,  it  cannot  be  matched  to 
damage  radii,  thus  reducing  the  mine's  effectiveness.  Consequently  some 
examination  of  the  effect  of  ferromagnetic  mine  cases  on  the  apparent 
sensitivity  of  magnetic  sensors  is  necessary. 


In  this  report  the  effects  on  sensitivity  to  an  external  field  of  a 
magnetic  sensor  placed  in  or  near  the  wall  of  a  ferromagnetic  case  are 
examined.  The  consequences  for  the  sensitivity  of  magnetic  mines  are 
discussed.  The  magnetic  field  problems  caused  by  ferromagnetic  mine  cases  are 
Investigated  using  analytic  and  numerical  techniques  in  two  dimensions 
(2D).  The  use  of  three-dimensional  (3D)  calculations,  while  possible,  would 
Involve  vastly  greater  computing  resources  and  the  development  of  codes  which 
may  not  prove  very  useful. 
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2. 


THEORY  OF  MAGNETOS TAT ICS 


Using  the  definitions  and  conventions  of  til,  the  basic  equations  of 
magnetostatics  are 
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which  reduce  to 


7.B  -  0 
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(2) 


in  the  absence  of  a  current  density  (J  -  o ) .  Equation  (2)  implies  that  the 
magnetic  field  intensity  (H)  is  derivable  from  a  scalar  potential 
♦  defined  by  H  •  -  7*.  Furthermore,  for  a  linear,  isotropic  material  the 
magnetic  Induction  B  -  where  »  is  the  magnetic  permeability  of  the 
material-  Then  7.B  «  o  implies  v.H-o  and  therefore  ♦"'satisfies 
Laplace's  equation  t 


7  ♦  -  o. 


(3) 


Laplace's  equation  is  commonly  encountered  in  electrostatics,  fluid  flow  and 
other  physical  problems.  In  2D,  Laplace's  equation  can  be  solved  analytically 
for  some  special  cases,  including  the  circular  tube  case,  but  in  general  the 
solution  is  obtained  numerically. 


Mo  problem  Involving  partial  differential  equations  is  completely 
specified  without  boundary  conditions.  The  problems  of  interest  here  Involve  a 
ferromagnetic  shell  (representing  a  mine  case)  in  an  otherwise  uniform 
magnetic  field  of  infinite  extent  (e.g.  the  problem  of  a  cylindrical  tube  of 
high  permeability;  Figure  1).  Two  types  of  boundary  conditions  arise 
naturally  from  3uch  a  problem.  Firstly,  the  field  far  from  the  anomaly  will* 
be  undistorted.  In  this  case  B  and  H  are  uniform  and  point  in  the  direction 
of  the  applied  field.  The  scalar  potential  *  and  the  field  B  are  given  by 


♦  -  -  B  rcose 
o 


(4) 


where  bq  is  the  magnitude  of  the  uniform  field,  r  and  a  define  the  location  of 
a  point  in  polar  coordinates,  and  x  is  a  unit  vector  parallel  to  the  uniform 
field.  The  other  boundary  condition  of  Interest  in  these  problems  is  at  a 
magnetic  interface  where 


i  irst  of  these  arises  from  the  conservation  of  magnetic  flux  across  the 
face  i.e.  there  are  no  sources  or  sinKs  of  magnetic  flux  (magnetic 
mo..opoles);  while  the  second  arises  from  the  continuity  of  the  tangential 
component  of  the  magnetic  field  Intensity  in  the  absence  of  a  surface  current. 


3.  LIMITATIONS  OF  THE  PHYSICAL  MODEL 


«* 

In  section  2,  an  approach  to  magnetostatics  problems  has 
been  developed  which  will  be  applied  to  specific  problems  in  sections  4,  5  and 
6.  Before  examining  these  specific  problems,  the  limitations  of  the  approach 
due  to  the  nature  of  magnetic  materials  will  be  examined. 


All  the  materials  under  consideration  are  assumed  to  have  linear  and 
isotropic  magnetic  properties.  Ferromagnetic  materials  exhibit  very  large 
magnetization  due  to  alignment  of  electron  spins  in  small  regions  Known  as 
domains,  which  in  turn  can  be  collectively  oriented,  -The  maximum  possible 
magnetization  is  obtained  when  all  available  electron  spins  are  aligned.  This 
is  Known  as  saturation  and  is  the  principal  cause  of  non-linear  behaviour  in 
magnetic  materials.  The  magnetization  of  a  given  sample  is  not  only  due  to 
the  external  field  but  also  to  the  material's  previous  history.  Thus 
ferromagnetic  materials  are  highly  non-linear  and  the  external  field  does  not 
completely  determine  their  magnetization.  Furthermore,  few  ferromagnetic 
materials  are  isotropic.  The  equation  B  -  *H  defines  >  rather  than  being  a 
material  constant  relating  B  and  H. 

The  assumptions  in  section  2  obviously  limit  the  usefulness  of 
performing  simple  magnetostatic  calculations.  However,  a  mine  case  in  the 
earth's  magnetic  field  is  very  unliKely  to  be  subjected  to  an  external 
magnetic  field  strong  enough  to  cause  the  case  material  to  approach  magnetic 
saturation.  The  normal  situation  of  interest  involves  a  small  transient 
change  to  a  larger,  effectively  constant  geomagnetic  field.  The  assumption  of 
linear  response  to  a  small  transient  will  be  reasonable  in  many  circumstances. 


As  already  discussed  permeability  is  not  a  well  defined  quantity  for 
ferromagnetic  materials.  The  Incremental  permeability  at  a  given  point  on  the 
hysteresis  curve  is  needed  for  a  given  problem.  In  general,  this  is  unKnown; 
however  usually  some  estimate  or  upper  bound  is  available. 


Given  these  apparent  limitations,  one  can  only  hope  to  obtain 
indicative  information  about  the  field  in  and  around  a  ferromagnetic  object 
without  extraordinary  information  about  the  magnetic  history  and  state  of  that 
object. 


4. 


ANALYTIC  SOLUTION  OK  CIRCULAR  TUBE  CASE 
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An  analytical  solution  of  the  circular  tube  case  (Figure  l)  is 
presented  permitting  a  test  of  numerical  solutions. 

The  solution  of  Laplace's  equation  for  a  cylindrical  problem  can  be 
solved  using  any  complete  set  of  functions.  Following  til,  for  a  cylindrical 
problem  the  natural  solution  for  the  potential  is  of  the  form 
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where  (a  )  and  (b  )  define  the  solution  in  a  region  and  the  P ^ ( cos^ )+^ye 
Legendre  polynomials .  Naturally,  as  r->0  the  terms  of  the  form  r 
diverge.  In  the  limit  r->-  the  the  non-vanishing  terms  must  tend  to  the 
solution  for  a  uniform  magnetic  field  i.e.  #  -  -  B  rcose. 

This  is  essentially  the  boundary  condition  far  from  the  anomaly  source  (4). 
Thu3  the  solution  must  be  of  the  form:- 
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where  (a  ),  if  ),  (7  )  and  (4  )  are  determined  by  the  boundary  conditions  at* 
the  magnetic  Interfaces  (5).  These  boundary  conditions  mean  that  only  terms 
with  i-l  can  be  non-zero.  Thus  we  can  drop  the  summations  over  t  and  the 
subscripts  1. 


From  (5)  the  boundary  conditions  at  the  magnetic  interfaces  (i.e.  at 
r-a  and  r-b)  give:- 
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which  can  be  solved  to  give 
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which  are  purely  functions  of  geometry,  f  and  the  external  uniform  field.  The 
magnetic  field  H  -  -  V#  is  readily  obtained  from  equations  (7)  and  (9)  as; 
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where  the  fields  have  been  transformed  from  cylindrical  polar  coordinates  into 
rectilinear  coordinates  with  x  -  rcost  and  y  «  rsin*. 

The  magnetic  lines  of  force  around  the  circular  tube  computed  from 
the  analytic  solution  are  shown  in  Figure  2  for  a-l,  b-2  and  0-1000.  The 
value  p-iooo  is  chosen  because  it  is  a  value  typical  of  iron-like  substances 
(21.  Figure  2  3hows  that  there  is  a  strong  concentration  of  flux  through  the 
ferromagnetic  material  ^th  the  highest  concentration  being  along  the  axis 
perpendicular  to  the  direction  of  the  external  field  (the  perpendicular 
axis).  The  magnetic  lines  of  force  are  almost  perpendicular  to  the  magnetic 
Interface  from  outside  at  r-b.  In  the  limit  0  -  they  would  be 
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perpendicular.  The  concentration  of  flux'through  the  ferromagnetic 
material  causes  increases  in  the  field  strength  outside  the  case  on  the 
magnetic  axis  and  a  decrease  along  the  perpendicular  axis.  These  distortions 
of  the  field  outside  the  case  constitute  the  magnetic  anomaly  and  have  short 
range,  dropping  off  as  1/r3  from  equations  ( x o ) .  inside  the  case  wall  is  a 
highly  shielded  region. 


5.  NUMERICAL  SOLUTIONS 


There  are  two  commonly  used  numerical  methods  for  the  solution  of 
second  order  partial  differential  equations;  the  finite  element  method  and  the 
boundary  element  method. 


In  view  of  the  lower  complexity  and  lesser  computational  power 
required,  the  boundary  element  method  is  used  to  examine  some  of  the  features 
of  the  magnetic  field  in  and  near  ferromagnetic  ob jeers.' 

if 

A  2D  boundary  element  code  written  in  FORTRAN  by  Chang  [3]  from  the 
University  of  Southampton  was  used.  The  code  was  developed  to  examine  the 
seepage  of  water  through  rocks  of  differing  permeabilities  under  a  given  head 
of  pressure.  This  problem  is  analogous  to  the  "flow"  of  magnetic  flux  through 
materials  of  differing  magnetic  permeability  under  given  initial  values  of 
*  at  some  boundaries.  The  only  changes  required  to  the  program  to  run  on  the 
MRL  VAX  11/780  were  file  input/output  statements.  In  order  to  get  best  use 
and  utility  from  the  code,  preprocessing  and  postprocessing  programs  have  been 
written  which  simplify  the  entry  of  data  to  the  boundary  element  program  and 
allow  various  types  of  data  output.  The  preprocessing  program  enabled  the 
rapid  generation  of  the  input  file,  which  contains  information  for  each 
boundary  element,  from  a  brief  specification  of  the  parameters  of  each 
boundary  line.  The  postprocessing  program  was  used  mainly  for  the  production 
of  diagrams  of  magnetic  lines  of  force.  Magnetic  lines  of  force  are  analogous 
to  streamlines  for  the  water  seepage  problem.  These  diagrams  provide  a  compact 
graphical  representation  of  the  general  characteristics  of  the  magnetic  field" 
in  a  region. 


6 .  RESULTS 


In  order  to  assess  the  value  of  numerical  solutions  to  the  2D 
magnetostatic  problems  under  study,  the  results  from  the  analytical  solution 
of  the  cylindrical  tube  problem  (section  4)  were  compared  with  those  from  the 
boundary  element  method  outlined  in  section  5. 

Numerical  solutions  were  obtained  using  45,89,178  and  338  boundary 
elements.  The  inner  radius  of  the  tube  was  1  unit,  the  outer  radius  2  units 
and  the  extremity  of  the  region  considered  was  10  units  except  when  338 
boundary  elements  were  used  when  this  was  increased  to  20  units.  The 
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arrangement  of  boundary  elements  for  the  numerical  solution  using  178  boundary 
elements  is  shown  in  Figure  3.  All  elements  were  made  to  have  approximately 
the  same  length. 

The  numerical  solutions  indicate  similar  gross  structure  to  the 
analytic  solution  which  can  be  seen  from  the  diagrams  showing  magnetic  lines 
of  force  (Figures  2£4). 

Figures  5-1 o  show  comparisons  of  the  analytic  solution  with  the 
various  numerical  solutions  using  differing  numbers  of  boundary  elements. 

The  fine  structure  in  the  analytic  solution,  S.g.  the  Increase  in  Hx  near  the 
tube  (Figure  6),  is  not  reproduced  by  the  numerical  solutions.  This  feature 
is  common  to  all  of  the  numerical  solutions.  Furthermore,  the  numerical 
solutions  do  not  asymptote  to  the  analytic  solution  with  increasing  numbers  of 
boundary  elements.  However  there  are  edge  or  boundary  effects  with  the 
boundary  element  method  which  are  related  to  the  size  of  the  elements. 
Increasing  the  number  of  boundary  elements  allows  solutions  to  be  obtained 
nearer  the  boundary  surfaces. 

The  implicit  averaging  in  the  boundary  element  method  seems  to  wash 
out  fine  structure.  This  does  not  poie  a  great  problem  considering  the 
inherent  limitations  of  the  physical  model  (section  3).  The  boundary  element 
method  provides  the  directions  of  streamlines  and  estimates  of  the  field 
strength  both  inside  the  material  and  in  the  tube  wall. 

Fortunately  the  poor  results  near  the  boundary  surfaces  are  not 
important  in  the  current  examination  of  mine  cases.  However  if  one  wished  to 
perform  calculations  on  a  comparatively  thin-walled  vessel  HJce  a  real  mine 
case,  the  boundary  effects  would  mean  that  many  elements  would  be  needed, 
greatly  Increasing  the  cost  of  the  computation. 

Figures  n-14  show  magnetic  lines  of  force  for  simple  2D 
ferromagnetic  objects  of  permeability  *  -  looo  in  an  otherwise  uniform 
magnetic  field.  These  objects  were  modelled  in  a  region  of  side  length  10 
units  with  boundary  elements  of  length  as  close  as  possible  to  those  in  the 
cylindrical  tube  example  using  178  boundary  elements  (Figure  3).  All  ► 
streamline  diagrams  are  at  the  same  scale. 

The  anomaly  properties  outside  the  object  are  generally  similar  to 
the  cylindrical  tube  example  with  an  increase  in  the  magnetic  field  strength 
along  the  axis  of  the  external  field  and  a  decrease  along  the  perpendicular 
axis.  In  the  case  wall  the  the  concentration  of  flux  again  has  similar 
properties  with  greatest  concentration  along  the  perpendicular  axis.  Inside 
the  objects  the  field  is  greatly  diminished  due  to  shielding. 

For  the  elliptical  cases  (Figures  14  and  is)  these  streamline 
diagrams  show  that  there  is  a  considerable  difference  in  sensitivity  for  mines 
aligned  with  their  major  and  minor  axes  parallel  to  the  external  uniform  field 
(ignoring  orientation  effects  e.g.  using  a  three-axis  magnetometer). 


CONCLUSIONS 


Although  these  calculations  were  performed  in  2D  for  simple  2D 
objects,  the  inferences  drawn  should  be  capable  of  extension  to  3D.  The 
diagrams  showing  magnetic  lines  of  force  clearly  indicate  that  the  apparent 
sensitivity  of  a  magnetic  sensor  will  be  greatly  decreased  if  it  is  placed 
inside  a  ferromagnetic  case.  If  the  sensor  is  placed  in  the  wall  of  a 
ferromagnetic  case  the  sensitivity  of  the  mine  will  depend  on  the  orientation 
of  the  case  relative  to  the  external  field. 

* 

variability  in  the  apparent  sensitivity  of  a  mine  sensor  makes  its 
use  much  more  difficult  in  that:- 
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The  sensor  response  to  a  given  target  cannot  be  accurately 
specified.  Consequently  valid  targets  may  be  missed  in  some 
orientations  and  damage  radii  cannot  be  properly  matched  to  target 
signal  with  a  consequent  reduction  in  the  eff ectlveness  of  one's 
mines . 


Difficulty  of  clearing  one's  odrn  minefields  is  increased  due  to  the 
need  to  use  worst  case  values  for  sensitivity  which  may  be 
unrealistic . 


Apart  from  these  operational  problems,  the  magnetic  anomaly  from  a 
ferromagnetic  case  provides  the  enemy  with  a  means  of  detecting  a  mine  that 
could  be  denied  him  if  a  non-magnetic  case  were  used. 

The  boundary  element  method  is  capable  of  being  used  to  determine 
relative  sensitivity  of  magnetic  sensors  Inside  mines  of  different  geometries 
given  that  other  parameters  (e.g.  orientation,  magnetic  history  of  material) 
are  unchanged.  A  knowledge  of  the  maximum  sensitivity  of  a  mine  is  necessary 
for  mine  clearance  operations.  An  estimate  of  the  maximum  sensitivity  of  a 
mine  may  be  obtained  from  a  boundary  element  calculation. 

The  actual  sensitivity  of  a  given  mine  cannot  be  predicted  if  the 
sensor  is  in  the  wall  of  a  ferromagnetic  case.  To  put  a  sensor  inside  the  case 
means  that  a  higher  sensitivity,  more  expensive  sensor  is  required  due  to  the 
shielding  effect  of  a  ferromagnetic  case.  The  use  of  a  ferromagnetic  mine 
case  means  that  variability  of  sensor  response  and  its  associated  problems 
are  accepted  by  the  designer  of  the  mine. 

Clearly,  in  the  design  of  a  mine  utilising  a  magnetic  sensor,  a  non¬ 
magnetic  case  and  a  three-axis  or  total  field  magnetometer  would  eliminate  all 
orientation  and  sensitivity  of  magnetic  sensor  difficulties. 
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FIGURE  1 


Cylindrical  tube  of  high  permeability  in  a  region  of  otherwise 
uniform  magnetic  field. 
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FIGURE 


Comparison  of  numerical  and  analytic  solutions  of  the 
cylindrical  tube  problem  along  the  line  x«y,  inside  and  in  the 
wall  of  the  tube. 
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FIGURE  10  Comparison  of  numerical  Qnd  analytic  solutions  of  the  cylindrical  tube  problem  along  the  axis  of  the  tube 
perpendicular  to  the  field,  outside  the  tube. 


FIGURE  14 


Lines  of  force  derived  from  the  184  b.e  numerical  solution  of 
the  elliptical  cross-section  tube  problem,  where,  the  major  axis 
of  the  tube  Is  aligned  with  the  field. 
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ABSTRACT 


In  order  to  assess  the  effect  of  a  magnetic  mine  case  on  a  magnetic 
sensor  in  or  near  the  case,  some  2D  magnetostatic  calculations  have  been 
performed  using  the  boundary  element  method.  For  the  case  of  a  cylindrical 
tube  the  analytic  solution  to  the  problem  is  presented  allowing  the  analytic 
and  numerical  solutions  to  be  compared.  Conclusions  are  drawn  about  the 
likely  effect  of  a  ferromagnetic  case  on  a  sensor  in  or  near  the  case. 
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